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Suppression of baryon number violation in electroweak collisions: Numerical results. 
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We present a semiclassical study of the suppression of topology changing, baryon number violating 
transitions induced by particle collisions in the electroweak theory. We find that below the sphaleron 
energy the suppression exponent is remarkably close to the analytic estimates based on a low energy 
expansion about the instanton. Above the sphaleron energy, the relevant semiclassical solutions 
have qualitatively different properties from those below the sphaleron: they correspond to jumps 
on top of the barrier. Yet these processes remain exponentially suppressed, and, furthermore, the 
tunneling exponent flattens out in energy. We also derive lower bounds on the tunneling exponent 
which show that baryon number violation remains exponentially suppressed up to very high energies 
of at least 30 sphaleron masses (250 TeV). 
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A long-standing problem in the electroweak theory 
is whether instanton-like processes occur at high rates 
in particle collisions near and above the sphaleron en- 
ergy. The energy of the sphaleron Q] represents the min- 
imum height of the barrier separating topologically dis- 
tinct vacua in a non-Abclian gauge-Higgs theory, thus it 
sets a non-perturbative energy scale at weak coupling. 
Instanton-like transitions between these vacua, which at 
low energies occur via tunneling and hence at exponen- 
tially small rates, are energetically allowed to proceed 
classically at energies above the sphaleron energy. The 
problem is whether or not these classical, and hence un- 
supprcssed transitions are allowed dynamically in colli- 
sions of highly energetic particles. This problem is par- 
ticularly interesting in the context of the electroweak the- 
ory, both because instanton-like transitions are accom- 
panied by non-conservation of baryon and lepton num- 
bers 0], and because the sphaleron energy is relatively 
low, E sph ~ 8TeV. 

As was first found in Rcfs. 0, 0, cross sections 
of collision-induced instanton processes increase rapidly 
with energy at E <C -Esph- Subsequently, it was 
shown 0, 0, S 13 that the total cross section has the 
exponential form 1 (for reviews see Refs. IE3EI) 

atot{E) ~ exp (-— F H G(E/E sph )\ , 
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where aw — ff 2 /47r is the small coupling constant (aw — 
1/30 in the electroweak theory). Perturbative calcula- 
tions about the instanton enable one to evaluate Fhg as 
a series in fractional powers of E/E sv ^ but the pertur- 
bative expansion becomes unreliable at E ~ -E S ph and at 
higher ener gies . Existing analytical estimates of Fhg at 
all energies |l2lll3j| are based on a number of assumptions 
which may or may not be valid. 

One way to understand instanton-like processes at high 
energies is to obtain numerically solutions to classical, 
real time field equations exhibiting appropriate topol- 
ogy and in this way explore the region of param- 
eter space where classical over-barrier transitions do oc- 
cur. Besides the total energy E, an important parameter 
is the number of incoming particles TV, which one cal- 
culates semiclassically for every solution. This approach 
enables one to find the approximate boundary of the clas- 
sicall y a llowed region in the (E, TV) plane; the analysis of 
Ref. [14J extends to E ~ 2£ , sp h and shows that even at 
the highest energy attained in this study the number of 
incoming particles is always large, TV ~ I /aw, which is 
very far from realistic collisions. 

In this paper we present the results of another com- 
putational approach, which is appropriate for analyzing 
the classically forbidden region in the (E, TV) plane. We 
study the four-dimensional SU (2) gauge theory with a 
Higgs doublet <f>, which corresponds to the bosonic sec- 
tor of the Electroweak Theory with 6w = and captures 
all relevant features of the Standard Model (to leading or- 
der, the effects of fermions on the dynamics of the gauge 
and Higgs field can be ignored The action of the 
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model is 
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In most of our calculations the Higgs self-coupling A was 
set equal to A = 0.125, which corresponds to Mr = Mw 
We found that the dependence of our results on the Higgs 
boson mass is very weak, so the specific choice of A does 
not affect our conclusions. 

Our starting point is the observation 0, 0, 0] that 
the inclusive probability of tunneling from a state with 
fixed energy E and fixed number of incoming particles N 
is calculable in a semiclassical way, provided that E = 
E I aw and N = N/aw, where aw is a small parameter 
and E and N are held fixed in the limit aw — > 0. This 
inclusive probability is defined as follows, 



a{E, N) 
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where S is the S'-matrix, Pe,n are projectors onto sub- 
spaces of fixed energy E and fixed number of particles 
N, and the states |i) and |/) are perturbative excita- 
tions about topologically distinct vacua. In the regime 
aw —* 0, with E and N held fixed, this probability can 
be calculated in the semiclassical approximation, leading 
to 



<r(E, N) ~ exp 
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where the exponent F(E, N) is obtained by solving a 
classical boundary value problem 0,0,0] about which 
we will have more to say later. 

Furthermore, it has been conjectured fl^. ll7Lll^ | that 
the exponent for the two-particle cross section is recov- 
ered in the limit of small number of incoming particles, 



Fhg(E) 



lim F(E, N) 



(2) 



This conjecture was checked in several orders of pertur- 
bation theory in E/E^ in gauge theory 0,^1 an d by 
comparison with the full quantum mechanical solution in 
a model with two degrees of freedom (2^, El ■ Hence, our 
strategy is to evaluate numerically F(E, N) in as large 
region of the (E, N) plane as possible, and then extrapo- 
late the results to N — > 0. In what follows we omit tilde 
over E and N to simplify notations. 

The boundary value problem for F(E,N) was derived 
elsewhere [lfjlllilGjL so we only present its formulation. 
Let cp(x,t) denote collectively all physical fields in the 
model. One introduces two auxiliary real parameters T 
and 9 and considers tp(x,t) as complex functions on the 
contour ABCD in the complex time plane shown in Fig.^ 
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FIG. 1: The contour in complex time used in the formulation 
of the boundary value problem . 



The parameter T determines the height of the contour 
(equal to T/2), while the role of 9 will be described later. 
The field <p should satisfy the field equations, 
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on the contour ABCD. In the infinite future (part D of 
the contour), the field should be real 

lnnp(x,Tf — > oo ) — > , 1m if (x,Tf — > oo) — > (3b) 

(for complex fields, such as <f> in QJ, this means that both 
(<f> + $ t )/2 and (§ — &)/2i must be real). The remaining 
boundary conditions are imposed in the infinite past, t = 
iT/2 + Ti, Ti — > — oo, part A of the contour. Since for 
Ti — > — oo the system reduces to a superposition of non- 
interacting waves about one of the gauge theory vacua 
(which we choose to be the trivial one for definiteness), 
the field ip linearizes 
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The boundary condition in the infinite past is then the 
u 9 boundary condition" 
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(3c) 



For 9 different from zero this equation implies that the 
fields themselves must be continued to complex values. 
For a complex field, like $ in (JJJ, its real and imaginary 
parts must be continued to complex values separately. 
Finally, there arc two more equations, 
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These equations indirectly fix the values of T and 9 for 
given energy and number of incoming particles. Note 
that they are in fact semiclassical expressions for E and 
N in terms of the frequency components of the incoming 
field. 
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Given a solution to the boundary value problem, the 
exponent for the inclusive transition probability is 



47T 

—F(E,N) 
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Prom Eq. (0J the variables (T, 0) appear to be Legcndrc 
conjugates to (E,N). This correspondence is strength- 
ened by the following relations. 
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These relations are useful as a cross check of the nu- 
merical procedure, and also as a mean of extrapolating 
F(E, N) to N = 0. 

This method of obtaining the exponent for tunneling 
probability was implemented in quantum mechanics of 
two degrees of freedom [53, El and in scalar theory 
exhibiting collision-induced false vacuum decay [2^]. It 
has been adapted to systems with gauge degrees of free- 
dom in Ref. [24j where preliminary study of the energy 
region below -E S ph was performed. 

Two remarks are in order. First, the boundary value 
problem (JSJ by itself does not guarantee that its solution 
interpolates between topologically distinct vacua. Ensur- 
ing that the solutions have correct topology is an inde- 
pendent and important part of the computational proce- 
dure. 

Second, we look for solutions to the boundary value 
problem which are spherically symmetric in space. 
Physically, since both instanton and sphaleron have this 
property, it is likely that the relevant solutions are also 
spherically symmetric. Technically, spherical symmetry 
reduces the number of equations considerably, so that the 
numerical analysis simplifies significantly. In the gauge 
Aq = 0, spherically symmetric configurations |25l | are 
parameterized by five two-dimensional fields a, a, (3, (i 
and v, 
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&(x, t) = [/x(r, t) + ii/(r, t)a ■ n]£ 



where £ is a unit two-column. The fields a, a, [3, /1, v are 
real in the original SU (2)-Higgs theory, but they become 
complexified due to the ^-boundary condition (|3c|l . 

We solved the boundary value problem J3J) numerically 
in the Aq = gauge on a grid of spatial size in radial 
direction R = 8/(\/2Mw) and number of spatial grid 
points N r = 90. The length of initial Minkowskian part 
of the contour AB was equal to 6/ (\2Mw)- The number 
of time grid points on this part was N t = 200 while on 
the Euclidean part BC it was equal to 150. The num- 
ber of time grid points on the part CD varied, with the 
maximum number of about 400. 
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FIG. 2: Lines of F(E,N) = const. The lines are labeled 
by the values of the suppression exponent —aw log a = 4ttF. 
The line labeled by (F = 0) is the boundary of the classically 
allowed region. The "fuzzy" line represents the approximate 
boundary of the classically allowed region found in the over- 
barrier calculations of Ref. [13 . 



The details of our numerical procedure are given else- 
where [2(| . Here we concentrate on our results. Clearly, 
only a part of (E, N) plane is accessible to numerical 
study: the difficulty of the calculation increases at higher 
energies and smaller number of particles, as the solutions 
get sharper and linearize slower at large negative times. 
The region of (E, N) plane covered in our study is shown 
in Fig. [3 where we present the results for F(E, N). Be- 
fore discussing the tunneling exponent F(E,N), let us 
comment on various types of solutions we have found. 
The upper left line is the line of periodic instantons. 
These are solutions to the boundary value problem with 
= [2H l2Sj , which correspond to transitions with the 
smallest tunneling exponent F for given energy. The line 
of periodic instantons ends at the sphaleron point 2 . The 
almost vertical line beginning at the sphaleron in Fig. |3 
separates two parts of the classically forbidden region in 
which the solutions have qualitatively different proper- 
ties. To the left of this line, the solutions are real on 
the Minkowskian part CD of the contour, and rapidly 
dissipate at large times forming spherical waves. This is 
illustrated in Fig. |3| where the field ( = (3 — ia is shown 
for a solution with relatively low energy 3 . On the other 
hand, in the right part of the forbidden region, the so- 
lutions are complex on the part CD of the contour, and 
obey the reality condition (|3b|) only asymptotically. Part 



The number of incoming particles for the sphaleron is obtained by 
infiiiitcsimally perturbing the (unstable) sphaleron solution alon g 
the negative mode, and integrating backwards in real time 1141 . 
Note that (3 = — 1, a = corresponds to the trivial gauge theory 
vacuum, while in the first topological vacuum £ winds around 
the unit circle in complex plane as r runs from to 00, sec Eg. 171 
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FIG. 3: The field ( for a solution with N = 1/aw and 
E — 3.35Mw /aw- The color tracks the phase of the field, 
whose behavior indeed shows the correct topology. The part 
corresponding to the Euclidean evolution is inclined for visu- 
alization purposes. 




FIG. 4: Field ( for a solution with N = 1/aw and E 
4.48Mw /a W - 



of the energy is emitted away in the form of spherical 
waves, but there remains a lump of energy near the ori- 
gin, as illustrated in Fig. 0] We have checked that this 
remaining lump is nothing but the sphaleron . 

The physical interpretation of the two types of solu- 
tions is as follows. At low energies, the system tunnels 
directly to the neighboring gauge theory vacuum plus 
linear excitations above it. At energies higher than the 
sphaleron energy (precisely, to the right of the almost 
vertical line in Fig. |2J) , the system ends up close to the 
sphaleron, with extra outgoing waves in the sphaleron 
background. In the latter case the system jumps on top 



4 The very fact that the solution is real only asymptotically in time 
is due to the existence of the unstable mode about the sphaleron. 



of the barrier. This process is not precisely what is usu- 
ally meant by tunneling; still it occurs with exponen- 
tially small probability which may be attributed to the 
rearrangement of the system from a collection of highly 
energetic incoming waves to the soft lump of the fields 
given by the sphaleron. The method of obtaining solu- 
tions of the second type was proposed in Ref. |22| . and 
we make use of this method in our work (see Ref. [26j for 
details). 

Let us now concentrate on the results for the tunnel- 
ing exponent F(E,N). Our data are in agreement with 
analytical results for F(E, N) in the low energy region, 
see Ref. [2^ for details. Another interesting comparison 
can be made with the results of Ref 01 > where a Monte- 
Carlo technique was used to find real-time overbarrier 
solutions close to the boundary of the classically allowed 
region. This technique produced an approximation (and, 
at the same time, an upper bound) for the boundary of 
the classically allowed region. It is seen that the results 
of Ref. are reasonably close to the boundary of the 
classically allowed region found in our calculations. 

To get insight into the suppression factor Fhg(E) for 
actual particle collisions, we have to extrapolate our data 
to N = 0, see Eq. We present here two types of 

extrapolation. The first one produces lower bounds on 
the suppression exponent Fhg(E) itself, while the sec- 
ond one gives an estimate for Fhg(E). While the latter 
extrapolation has stronger predictive power at relatively 
low energies, E < 2E sp h, the former extends to much 
higher energies, so the two are complementary. 

We begin with the lower bounds on Fhg(E). One 
way to obtain a lower bound is to make use of Eq. JSJl, 
together with the fact that 9 increases as N gets smaller. 
Hence, a lower bound on Fhg(E) is obtained by simply 
continuing F(E, N) with a linear function of N for each 
energy. This bound is shown in Figs. 03 El dashed line. It 
indicates that up to the energy %Mw/otw — 20TeV the 
suppression is still high: the suppression factor is smaller 
than e" 60 - 1(T 26 for a w ~ 1/30. 

Another lower bound, the best we can obtain at very 
high energies, is constructed by exploiting the observa- 
tion that the lines of constant F in (E, N) plane have 
positive curvature (see Fig. |2J). So, the lower bound is 
obtained by extrapolating these lines linearly to N = 0. 
This bound is displayed in Fig. dashed-dotted line. 
One can see that exponential suppression continues up 
to the energy of at least 250 TeV. 

Let us now come to the second type of extrapolation 
which wc make to estimate Fhg(E) itself. We find it 
appropriate to use Eq. ©. The point is that the func- 
tion T(N) at fixed energy is approximately linear in N. 
This property has been shown analytically for low ener- 
gies |2!j , while for all energies it follows from our numeri- 
cal data. [This is in contrast to the behavior of F(E, N): 
the analytical results at low energy show that for fixed 
E, this function behaves as NlogN as N — > pflj.] 
We thus extrapolate T(N, E) linearly to N = along 
E = const, and then integrate Eq. JSJ at N = to ob- 
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FIG. 5: Lower bounds on the suppression exponent for two- 
particle collisions, dashed and dashed-dotted lines. The dot- 
ted line is the estimate of Refs. [l2l H^[. 
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FIG. 6: Estimate of the suppression exponent for two-particle 
collisions Fhg(E) (solid line), lower bound on Fhg(E) 
(dashed line), low energy analytic prediction (rare dot- 
ted line) and analytic estimate of Refs. |l2l Il3jl (dotted line) . 



tain the suppression exponent Fhg(E) for two-particle 
collisions. This estimate is shown in Fig. |SJ solid line. 
It is instructive to compare it to the one loop analytic 
result [3(], HOI HE 0- which gives three terms in the 
low-energy expansion, 



^F(E) = ±L 
aw otw 



1 ~8\~Eh) + 16\E^ 



(8) 



where Eq = \/Q-KMw/aw- We see that our numerical 
data are (somewhat unexpectedly) very close to the one 
loop result (JSJ) up to the sphaleron energy. In this en- 
ergy region, they are consistent also with the analytic 
estimate of Refs. ^E 113 • O n the other hand, the behav- 
ior of Fhg{E) changes dramatically at E > _E sp h. We 
attribute this to the change in the tunneling behavior — 
at E > -Esph the system tunnels "on top of the bar- 
rier" . Our numerical data show that the suppression ex- 
ponent Fhg(E) flattens out, and topology changing pro- 
cesses are in fact much heavier suppressed at E > -E S ph 
as compared to the estimate (JSJl and the estimate of 
Refs. |H|l3. 

Thus, our numerical results, albeit covering a limited 
range of energies and initial particle numbers, enable us 
to obtain both a lower bound for and an actual estimate 
of the suppression exponent for the topology changing 
two-particle cross-section in the electroweak theory well 
above the sphaleron energy This cross section remains 
exponentially suppressed up to very high energies of at 
least 250 TeV. In fact, the energy, if any, at which the 
exponential suppression disappears, is most likely much 
higher, as suggested by comparison of our lower bound 
and actual estimate at energies exceeding significantly 
i^ph, see Fig. El 
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